Leading twist distribution amplitudes of P-wave nonrelativistic mesons. 
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This paper is devoted to the study of the leading twist distribution amplitudes of P-wave non- 
relativistic mesons. It is shown that at the leading order approximation in relative velocity of 
quark-antiquark pair inside the mesons these distribution amplitudes can be expressed through 
one universal function. As an example, the distribution amplitudes of P-wave charmonia mesons 
are considered. Within QCD sum rules the model for the universal function of P-wave charmonia 
mesons is built. In addition, it is found the relations between the moments of the universal function 
and the nonrelativistic QCD matrix elements that control relativistic corrections to any amplitude 
involving P-wave charmonia. Our calculation shows that characteristic size of these corrections is 

0"N • of order of ~ 30%. 
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?H 1 I. INTRODUCTION 



Hard exclusive processes are very interesting both from theoretical and experimental points of view. Commonly, 
theoretical approach to the description of such processes is based on the factorization theorem [J Q • Within this 
| theorem the amplitude of hard exclusive process can be separated into two parts. The first part is partons production 
at very small distances, which can be treated within perturbative QCD. The second part is hardronization of the 
partons at larger distances. This part contains information about nonperturbative dynamic of strong interaction. For 
hard exclusive processes it can be parameterized by process independent distribution amplitudes (DA) , which can be 
considered as hadrons' wave functions at light like separation between the partons in the hadron. It should stressed 
Ph' that DAs are very important for the calculation of the amplitude of any hard exclusive process. 

Recently, the leading twist DAs of 5-wave nonrelativistic mesons have become the object of intensive study [H, 0, 
H, H0, Knowledge about these DAs allowed one to build some models for S'-wave charmonia DAs, that can 

be used in practical calculations. In this paper general properties of the leading twist DAs of P-wave nonrelativistc 
mesons will be studied. The results of this study will be used to build the model for the DAs of P-wave charmonia 
mesons, that can be used in calculations. 

This paper is organized as follows. In the next section all definitions of the DAs of P-wave mesons will be given. 
These DAs will be studied in section III at the leading order approximation in relative velocity of quark-antiquark 
■ pair inside P-wave meson. In section IV QCD sum rules will be applied to the calculation of the moments of the 
P-wave charmonia DAs. Using the results of this study the model of the P-wave charmonia DAs will be build in 
section V. In the last section the results of this paper will be summarized. 
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II. DEFINITIONS OF THE DISTRIBUTION AMPLITUDES. 

In this section the definitions of the leading twist distribution amplitudes (DA) of P-wave nonrelativistic mesons 
will be given. In the conventional quark model nonrelativistic mesons are quark-antiquark((5Q) bound states. In 
these mesons the quark-antiquark pair can be in the spin singlet or spin triplet states. Since orbital momentum of the 
quark-antiquark pair is unity one can conclude that there are four P-wave mesons: Xo( 3 -Po)7 Xi( 3 -Pl)j X2( 3 -P2), h( Pi). 
The leading twist DAs of these mesons can be denned as follows, 
for the xo meson: 

(Xo(p)\Q(zh»[z,-z}Q(-z)\Q) = f xoPfi f dte^<f> X0 {tri, (1) 
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for the xi meson: 

(Xi(p,ex= )\Q(zh ll j 5 [z,-z}Q(-z)\0) = f^ Plt J d{ p), 

(Xi(p,e x =±i)\Q(z)a flv [z,-z}Q(-z)\0) = /^e^eV J ' d£ e^<f> Xl (£ M ), (2) 

for the /i meson: 

(h(p 7 e x=±1 )\Q(z)<j^{z,-z}Q(-z)\0) = f' h e ltuafl e a pP J d£ <f>' h fo > ( 3 ) 

for the X2 meson: 

(X2(P, e x=0 )\Q(z) 7ll [z, - z ]Q(-z)\0) = f X2 p„ f d£ e^4> X2 (£, //), 



(Xa(p,eA=±i)|0(2)^[z,-«]Q(-«)|0) = f X2 M X2 ( PfiPv - p uPfi ) [ dte^^&fx), p M = (4) 

where the following designations are used: 3Gi,X2 are the fractions of momentum of meson carried by quark and 
antiquark correspondingly, £ = X\ — x%, p, e are the momentum and polarizations of P-wave mesons. For the mesons 
Xi , h the polarization e is described by the four vector e M , for the X2 meson the polarization e is described by the 
tensor e M „. The factor [z, —z], that makes matrix elements (H])-® gauge invariant, is defined as 

[z,-z}= Pexp[ig f dx^A^x)}. (5) 

J — z 

In applications it is useful to rewrite the four-vector p M in the following way. Evidently, one can write the polarization 
of the X2 meson in terms of the polarization of two vector mesons. Thus for the transverse polarization of the meson xi 
one has £ x =±i ~ ( e A=±i ' e A=o + e A=o ' e x=±i)/^ ^t,v e>lv — !)■ If we further contract the polarization tensor e M „ with 
lightlike four-vector z, to the leading twist accuracy we will get z v = e x _ ±1 {pz) / '(y/2M X2 ) or p M = e^ =±1 /(v / 2M X2 ). 
This form of the vector p can be used in the calculation with the leading twist accuracy. It should be noted that the 
states of the X2 meson with the polarizations A = ±2 give contribution only to higher twist DAs. Since, this paper is 
devoted to the study of the leading twist DAs we don't consider these states. 

The functions without primes <p x o(£), <f> x i(0> 4>h(Q, 0x2 0x2 (£) are £ °dd and they are normalized as 

J dZZ<Kt) = l. (6) 

The functions with primes 4>' x i{£), <fi'h(£) are £ even and they are normalized as 

1 d£ ^(0 = 1. (7) 
In this paper all DAs will be parameterized by their moments 

<r> = / # em- («) 



Evidently, for the DAs with primes all odd moments are zero. For the DAs without primes all even moments are zero. 
To separate the moments of different DAs of the %2 meson, below we are going to use the following designations: (£") 
for the moments of the function X2 (£), for the moments of the function <fi X2 (£). 

The distribution amplitudes 0(£) and the constants / that parameterizes corresponding currents (HJ)— (I3J) are scale 
dependent objects. For applications it is useful to write how they depend on the scale |2j. To do this we expand DA 
in the series 



r oo 



(9) 
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where C^J 2 {C) are Gegenbauer polynomials. At the leading logarithmic accuracy the coefficients a n are renormalized 
multiplicatively 

a " (Ai)= (^^) 7 " /6 ° a ™' T(M0) ' (10) 
where 7„ are the anomalous dimensions. For the current Q<J U u[z 1 —z]Q the anomalous dimcntions are 

In = Cf[l + AY^-y 6 = ll--n fl , C f = -, (11) 

for the other currents 

(12) 

It is clear that the DAs without primes contain only n-odd terms in series ©. DAs with primes contain n-even terms 
in the series. 

The constants defined in equations (jT])- (|4]) are multiplicatively renormalizable. Using formulas (|9j)- (fT2|) one can 
determine the evolution of these constants 




For the constants f Xo , fh, f X2 the anomalous dimentions 7 is equal to 8/3C/, for the constants f Xl ,f X2 7 = 3C/, 
for the constant f' h 7 = Cf, for the constant /' 7 = 0. It should be noted that if the anomalous dimensions 7 of 
the constants / are factored from sum ([5]) , the anomalous dimensions of the remaining terms equal to the difference 

In ~ 7- 



III. THE DISTRIBUTION AMPLITUDES AT THE LEADING ORDER APPROXIMATION IN 

RELATIVE VELOCITY. 

In this section the DAs under study will be considered at the leading order approximation in relative velocity of 
quark-antiquark pairs inside the mesons. First let us consider the DA of the xo meson. The moments of this DA can 
be represented as follows 

f Xo (pz) n+1 < C > X0 = ( X o\Q2(-izD) n Q\0). (14) 

To get the expressions for the constant f Xo and the moment < £" > Xn one needs to calculate the matrix element in 
the right hand side. At the leading order approximation this calculation can be done using projector [TTL [l2| 

0(p)Q(p)- / <k^j= A 1 -(p-m Q )r(P + 2E)(p + m Q ), (15) 

J x/SmQ AV2E(E + m Q ) 

where P, q are the total and relative momentum of the QQ pair, mQ is the mass of the quark Q,p = P/2+q,p = P/2—q, 
E 2 = P 2 /4 = rriQ—q 2 . The matrix r = 75, eg for the spin singlet and spin triplet quark-antiquark pair correspondingly 
where es is the spin polarization of this pair. The scalar products P ■ es = 0, P ■ q = 0. In the center mass frame the 
dq is reduced to the <i 3 q/(27r) 3 and the function (p(—q 2 ) is reduced to the usual nonrelativistic wave function 0(q 2 ). 
For the xo meson the wave function y(q 2 ) can be written in the form 

^(q 2 ) = (16) 

At the leading order approximation in relative velocity the function ^>(q) is universal function for all P-wave mesons. 
It is normalized as 

d 3 a 

%VWq)| 2 =^. (17) 



(27T) 
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With this normalization of the function ^>(q), the function y>(q) is normalized as 

The same normalization condition will be used for the wave functions of all mesons under consideration. Using 
equations (|T4"J) . (fT5|) and (fT6|) one gets the result 



V m<3 

where M Xo is the mass of the \o meson, A n equals to 



fxo < C +1 > xo = -2-+V- ^4,, (i9) 

w mQM£o n + 3 



' r( ' 1 'q|>(q). (20) 



(27T) 



It should be noted that because of the Coulombic part of the nonrelativistic QCD potential the right hand side of 
equation Q19p is ultraviolet divergent (l3l . [lij . The moments of the DA in the left hand side are QCD operators (see 
equation ifTS]) ). In full QCD the Coulombic part of the nonrelativistic potential corresponds to the rescattering of 
the quark-antiquark pair of the QCD operators which is also ultraviolet divergent 0. To control the divergences 
in the right and left hand sides of equation (I19p it is assumed that both sides are regularized within dimensional 
regularization. 

It is interesting to note that relation (jT9|) is closely connected with Brodsky-Huang-Lepage (BHL) [l5| procedure. 
This fact can be seen as follows. Let us rewrite this relation as follows: 

/«*»*»(£) ~ ~ ^ h 2 ^^ ~ JH^r^- (21) 

Note that q z ip{q) is the L z = component of the wave function y(q 2 ) which is the only component important for the 
leading twist DA. So, the last relation can be written as follows 

/ ^ n+ Vxo(0 ~ J ^V^=o(q 2 ). (22) 

Further, we change the variables in the right side of this equation 



„ M 2 + ql 

q±^q±, Qz^ZMo, M 2 = * X . (23) 

Note also that £ <C l,q_i_ Mq and at the leading order approximation in relative velocity of the quark-antiquark 
pair in the meson relation (|2"2")) can be written as follows 

/f f M 2 £ 2 + a 2 \ 

dHC +1 x J d 2 q ± <p La=fS f \_ e ) ■ ( 24 ) 



So, at this level of accuracy the DA is just 



M 2 ( 2 + q 2 



, Xo (0 ~ / d V <p Lx=0 ( f ^ ) , (25) 



what coincides with BHL procedure. 

It is not difficult to get the relations for the other DAs and mesons 

Xi meson: 0(q ) = a — ^(q) 

ft, meson: </>(q 2 ) = 6 ^(q.), 
v2 

^2 



X2 meson: 0(q ) = e. iJ e s 1 q J V'(q), (26) 
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where e is the polarization three vector of the \i an d h mesons, is the polarization tensor of the \2 meson. Using 
these expressions one can get 



fx, <C +1 > Xl =-2" +1 



3 A n+2 1 
m Q M™ +1 n + 3' 



m Q MZ+ 1 (n + l)(n + 3)' 
1 




m Q M™ +1 (n + l)(n + 3)' 



f ^ tn+l v, _ _on+l / 4 A„+2 1 



^<^ +i >-=- 2n+ V4^^' (2?) 

The moments with primes correspond to the DAs with primes. As it was noted above, since all the DAs without 
primes are £-odd only odd moments of these DAs are nonzero. Similarly, for the DAs with primes only even moments 
different from zero. For these reason we took (n + l)th moments of DAs without primes and nth moments with 
primes, where n is assumed to be even. It should be noted that since the function ^(q) is universal for the P-wave 
mesons, the A n is the same for all mesons under study. From relation P7|) one can conclude that all the constants fi 
can be expressed through the only constant which will be designated as F 



F = V3/ Xo = V2f xt =i^ = fh = if' h = = L ■ (28) 



From equations (|27p one can also find that at the leading order approximation the nth moments of all functions 
without primes coincide. This means that these DAs are equal to one universal DA which will be designated as 

*(0 = = <M0 = MO - <M0 = <M£). (29) 

The same is true for all DAs with primes, which are equal to one function which will be designated below as 

*(0 = <x (0 = ^(0- (so) 

In addition, one can relate the moments of the $(£) to the moments of the as follows 

< c >*= <r l l 1 > " - ( 31 ) 

n + 1 

It should be noted that the constants and DAs in relations (f!?5|) -(j3"T j) depend on scale in a different way. This means 
that relations (|28[) - (f5T|) . which are valid at not to large scale, will be violated at sufficiently large scale. 

Recursive relation (f3Tj) determines the function through the function $(£). One can guess the solution of this 
relation: 

(32) 



To prove that (|32|) is the solution of relations (|3T|) one should put this function to the definition of the n-th moment 
([5]) and integrate the resulting expression by parts. It is seen from equations \29fl , and that all DAs of the 
P-wave mesons are defined through the universal function 4 r (£). It should be noted that this fact results from the 
nonrelativistic spin-symmetry, which holds at leading order in the heavy-quark velocity. Below equations f|29[) and 
(f32ll will be used to build the models for the function <&(£) and \I r (£) of P-wave charmonia. 
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At the end of this section it is interesting to discuss the question about relativistic corrections to the matrix elements 
involving P-wave quarkonia. If we ignore the contribution coming from the higher Fock states, relativistic corrections 
to the matrix involving, for instance, the xo meson is given by the matrix element 

1 ( Xo | X +(aD)DV|0) 
( )P ~K>) 2 <XO|X + (^D)V>|0) » ^ 

where ip and x + are Pauli spinor fields that annihilate a quark and an antiquark respectively, a are Pauli matrixes, 
m*Q is the quark pole mass. Using equation (|19[) it is not difficult to obtain general formula that connects the moment 
with the matrix element (v n )p at the leading order approximation in relative velocity 

Ai+2 , ^, n +2\ _ n + ^ , cn +l 



0(v n+2 ) = -4— (r + ) + 0(v n+2 ). (34) 



Although the derivation was done for the \o meson, at the leading order approximation in relative velocity the matrix 
element (v n ) p is universal for all P-wave mesons. 



IV. CHARMONIA DISTRIBUTION AMPLITUDES WITHIN QCD SUM RULES. 

In this section QCD sum rules approach [l6|, QjJ will be applied to the calculation of the moments of the P-wave 
charmonia DAs. The problem which can be met if one tries to apply QCD sum rules to the calculation of these 
moments is that it is not possible to write QCD sum rules for one DA. Commonly, the contribution of different 
mesons and different DAs mix in QCD sum rules, what does not allow us to calculate the constants and the moments 
of DAs (|T|)-(j4j) separately. Only some combinations of different constants and moments can be extracted from QCD 

sum rules. For instance, two point QCD sum rules for the currents Q(x)j IM (zD)Q(x) ■ Q(0)j tJ/ (zD) n+1 Q(0) can be 
written as follows 

/L^Vo Zxc2^" +1 )xc2 1 f° Im n pert (5,n) (m) 2 

Kco + Q 2 ) m+1 Ke 2 + Q 2 ) m+1 ~ * hral (S + Q 2 ) m+1 U »* e * tW ' 

The expressions for the Im n pe rt(s, n), n^™] rt (<2 2 , n) will be given below (equations ([37]) . ([38]) ). It is seen from this 
example that it is not possible to extract the constants / Xc0 , /x„2 or the moments (£™ +1 ) Xc0 , (£™ +1 )x<:2 f rorn (EH) 
separately. Evidently, this strongly restricts the accuracy of the calculation. The only QCD sum rules which are free 

from this problem are two point sum rules for the currents Q{x)^^ l 'y^{zD)Q{x) ■ Q{Q)^^ 5 (zD) n+1 Q{Q) 

fhM n+1 )h c 1 f s ° Im n pcrt (s,n) (m) 2 

dS f.. fflm +l + n npcrt(Q ,»), (36) 



K e + Q 2 )™+! 7T J 4m 2 [S + Q 2 )™+! 

where Im n port (s, n), ilj™^ (Q 2 , n) can be written as 

n^ rt (g 2 ,n) = n^(g 2 ,n) + n^(Q 2 ,n) + n^(Q 2 ,n), (38) 

#»)rr» ^ Mluxn f 1 M , n(n+l)^ ri ^ (l-£ 2 r +2 



^\Q 2 ,n) = <^>„(„+i) (m +i) f\dte- fA , , _ _, m+2 . 



I I L „ \ I I L KJ lib \ I U L L I J. L/L I „,. . o . 

6^ V-i V ; (4m2 + Q2(l-^))"^' 
iH/n2 ^ _ ( a « G2 )„^ i iu m i -n [ cn (i 



(4m2 + Q 2 (l-e 2 )) 

In the calculation we take Q 2 = 4m 2 |18j. In the numerical analysis of QCD sum rules the values of parameters 
and (a s G 2 /n) will be taken from paper [18l |: 



m c = 1.24 ± 0.02 GeV, (— G 2 ) = 0.012 ± 30% GeV 4 . (39) 

7T 



7 



First sum rules (136]) will be applied to the calculation of the constant /jj . It is not difficult to express the constant 
f? from equation (J36J) at n = as a function of m. For too small values of m (m < mi) there is large contributions 
from higher resonances and continuum which spoil sum rules (|36|) . Although for m m\ these contributions are 
strongly suppressed, it is not possible to apply sum rules for too large m (m > 7712) since the contribution arising from 
higher dimensional vacuum condensates rapidly grows with m what invalidates our approximation. If ra\ < rri2 there 
is some region of applicability of sum rules (|36|) [mi, 7712] where the resonance and the higher dimensional vacuum 
condensates contributions are not too large. Within this region as a function of m varies slowly and one can 
determine the value of this constant. The value of the continuum threshold s$ must be taken so that to appear 
stability region [7711,7712]. Our calculation shows that for the central values of parameters (l39l) there exists stability 
region for sq > (4.3 GeV) 2 . If the value of the continuum threshold sq is varied in the region so £ (4.3 2 ,oo) GeV 2 , 
the value of the constant /| can be written as f 2 c = (0.037 ± 0.005) GeV 2 . In addition to the uncertainty due to the 
variation of the value of sq, there are uncertainties due to the variation of the values of the m c (which is ±0.004) and 
(a s /irG 2 ) (which is ±0.001). The last source of uncertainty is the radiative corrections to the perturbative density 
ImIIp Cr t(s, n), which will be estimated as a s (m c )/ir ~ 13%. Adding these uncertainties in quadrature, one gets 

/ 2 c (M~ m c ) = (0.037 ±0.007)GeV 2 . (40) 

As it was noted above, the value of the constant f% is scale dependent quantity. The characteristic scale of QCD sum 
rules is ~ m c . This means that the value of the constant f 2 is determined at the scale ~ m c , as it is shown in (|40|) . 

Next let us consider the moments (£ 3 )/i c , (£, 5 )h c , (£ 7 )/i c - However, instead of considering QCD sum rules for n = 2, 4, 6 
we will consider the ratios of sum rules at n — 2, 4, 6 and the sum rules at n = 0. Such approach improves the accuracy 
of the calculation (see paper [5j for details). The analisys similar to that for the constant gives 

(f) = 0.18 ±0.03, 
(£ 5 ) = 0.050 ±0.010, 

(£ 7 ) = 0.017 ±0.004. (41) 

It should be noted that these moments are defined at the scale ~ m c . It should be also noted that the values 
of the moments (|41 [) are in good agreement with potential model estimation (see paper Q). For instance, within 
Buchmuller-Tye potential model pj| (£ 3 ) = 0.18, (£ 5 ) = 0.047, (£ 7 ) = 0.016; within Cornell potential model [13] 
(£ 3 ) = 0.16, (£ 5 ) = 0.040, (£ 7 ) = 0.013. 

Using values (|4ip one can find the matrix elements that control relativisitic corrections to any process with P-wave 
charmonia in the initial or final state. The relationships between these matrix elements and the moments are given 
in equation (|34j) . These relations are valid up to the higher order relativistic corrections, which can be estimated as 
(v 2 ). Taking into account this additional source of uncertainty one gets 

(v 2 )p = 1(f) p = 0.30 ±0.10, 
(t/) p = I(e 5 ) P - 0.12 ±0.04, 

(v e ) P = ^(£ 7 ) P = 0.051 ±0.018. (42) 
In the next sections results ([4"T]) will be used to build the model for the DAs $(£) and '!'(£). 



V. MODEL FOR CHARMONIA DISTRIBUTION AMPLITUDES. 

To build the model of the function < I , (£,/7 ~ m c ) we use Borel version [l6|, [ljj of sum rules (|3"6"|) but without 
continuum contribution and power corrections 

/ina-< /M! = ^ Si r +2 |(1 - e 2 )exp(-^ T ^) . (43) 
Evidently, within this approximation the function Q(£,fi ~ m c ) can be written in the form 

m fi ~ m c ) = c(/? P )(l - £ 2 ) C ex P (-J^j , (44) 
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where c(/3p) is a normalization constant and (3p is some constant. We propose function (|44|) as the model for DAs 
^(Cj A 1 ~ m c)- To fix the constant (3p the value of the moment (£ 2 ) (|4"Tj) will be used. Thus we get (3p — 3.4to'g- The 
constant c(/3) can be determined from normalization condition (|6|). The moments of the function l[44|) are 

(£ 3 ) = 0.18 ± 0.03, (45) 
(£ 5 ) = 0.047 ±0.014, 
(£ 7 ) = 0.015 ±0.006. 

Using the model for the function /i ~ m c ) and equation (13211 one can get the model of the DA fi ~ m c ) 

p ~ m c ) = - ^ dt3(t, /x ~ m c ) = ^M(l - C 2 ) 2 ^3 (3^2) ' ( 46 ) 

where the function E^{z) is the exponential integral function 

/°° e~ zt 
dt—. (47) 

To determine DAs CO)-© at a scale different from m c one should use relations (J23), (f3T)]) and than evolution equations 
for the corresponding DAs. 



VI. CONCLUSION. 



In this paper we have considered the leading twist distribution amplitudes (DA) of P-wave nonrelativistic mesons. 
At the leading order approximation in relative velocity of quark-antiquark pair inside the mesons these functions can 
be expressed through one universal DA. We have derived the relations between the DAs of P-wave nonrelativistic 
mesons and the universal DA. 

As an example, we have considered the DAs of the P-wave charmonia mesons. Within QCD sum rules we found 
the moments of the leading twist DA of h c meson, what allowed us to build the model of the universal DA for the 
P-wave charmonia mesons. 

In addition, we have found the relations between the moments and the nonrelativistic QCD matrix elements that 
control relativistic corrections to any amplitude involving P-wave charmonia. The calculation shows that characteristic 
size of these corrections is ~ 30%. 
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